Enhanced magnetocaloric effect in frustrated magnets 
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The magnetothermodyriamics of strongly frustrated classical Heisenberg antiferromagnets on 
kagome, garnet, and pyrochlore lattices is examined. The field induced adiabatic temperature change 
(8T/dH)s is significantly larger for such systems compared to ordinary non-frustrated magnets and 
also exceeds the cooling rate of an ideal paramagnet in a wide range of fields. An enhancement of 
the magnetocaloric effect is related to presence of a macroscopic number of soft modes in frustrated 
magnets below the saturation field. Theoretical predictions are confirmed with extensive Monte 
Carlo simulations. 
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I. INTRODUCTION 

The magnetocaloric effect consists in heating or cool- 
ing of matter in response to variations of external mag- 
netic field. The technique of adiabatic demagnetization 
based on the magnetocaloric effect has been successfully 
utilized over the years to reach temperatures in a sub- 
Kelvin range.El Because of its technological simplicity and 
independence of gravity the magnetic cooling steadily at- 
tracts attention as a rival to the .dilution refrigerators, es- 
pecially in satellite applications. aWa Recent experimental 
studies of the magnetocaloric effect in specially designed 
materials have revealed the potential applicability of the 
magnetic coolina-|technique for room temperature refrig- 
eration as well.uEnD 

Paramagnetic salts, which are standard refrigerant ma- 
terials for the low-temperature magnetic cooling, contain 
dilute systems of magnetic dipoles. The entropy of an 
ideal paramagnet depends on magnetic field H and tem- 
perature T via H/T. Therefore, during adiabatic de- 
magnetization, temperature of a paramagnet decreases 
linearly with a field: T/H = const. The higher the den- 
sity of the magnetic moments and their spin number, the 
greater the cooling power of a refrigerant is. With in- 
creased density of spins, however, the strength of residual 
interactions between magnetic dipoles grows and leads 
to an ordering or a spin glass transition. The transition 
temperature T c , which is roughly given by a strength 
of magnetic interactions, limits the lowest temperatures 
achievable with paramagnetic salts. Even above T c the 
cooling becomes less efficient due to a reduced entropy, 
and the adiabatic temperature change is slower than in 
an ideal paramagnet. The cooling rate (dT/dH)s of a 
salt does not, therefore, exceed the corresponding rate of 
a paramagnet: (8T/dH) p s aia = T/H. 

The experimental and theoretical studies over the past 
decade have established a new distinct class of magnetic, 
materials called geometrically frustrated magnets JjuEll 
Magnetic ions in these systems form special types of 
lattices, Fig. [|. Despite interaction between neighbor- 
ing spins, strongly frustrated magnets remain disordered 
and possess finite entropy at temperatures well below the 



Curie- Weiss constant. This property suggests frustrated 
magnets as prospective candidates for use in the adia- 
batic demagnetization refrigerators. In fact, gadolinium 
gallium garnet Gd3Ga50i2 (GGGAhas long been known 
as a suitable refrigerant material.til Recently, GGG with 
its unique phase diagram in magnetic field attracted a lot 
of interest from the point of view of geometric magnetic 
frustration on a garnet lattice.E3 In the present work we 
demonstrate a close relation between an enhanced mag- 
netocaloric effect and strong geometric frustration. In 
particular, (dT/dH)$ in a frustrated magnet can exceed 
the cooling rate of an ideal paramagnet by more than an 
order of magnitude. 

Specifically, we study an unusual magnetothermody- 
namics of classical Heisenberg antiferromagnets on sev- 
eral frustrated lattices. The choice of classical models is 
motivated by large values of spins of the magnetic ions 
in frustrated magnetj&nCr 3+ (S = 3/2), Fe 3 + (S = 5/2), 
and Gd 3+ (S = 7/2). ffl The studied frustrated structures 
include two- (2D) and three-dimensional (3D) networks 
of corner-sharing triangles known as a kagome, Fig. 0a, 
and a garnet lattice, Fig. [j]b, and a 3D network of corner- 
sharing tetrahedra of a pyrochlore lattice, Fig. 0c. In all 
three cases the magnetic Hamiltonian includes exchange 
interaction of strength J between the nearest-neighbor 
classical spins of unit length and the Zeeman energy: 
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The main distinction of strongly frustrated magnets is a 
large, macroscopic degeneracy of their classical ground 
states in zero magnetic field. It arises due to a so called 
under-constraint: condition of the minimal energy for 
every elementary frustrated unit, triangle or tetrahe- 
dron, isjiot sufficient to fix all microscopic degrees of 
freedomJ§E£l In applied magnetic field the macroscopic 
degeneracy persists up to the saturation field H sat , which 
is -ff S at = 6 J for kagome and garnet antiferropiaffinets and 
-Hsat = 8 J for a pyrochlore antiferromagnetiijlia At zero 
temperature frustrated magnets loose their degeneracy 
and become completely polarized above H sa ±- Transfor- 
mation from the unique ground state at H > -ff sa t to 
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FIG. 1: Magnetic systems with strong frustration: (a) kagome lattice, (b) garnet lattice, and (c) pyrochlore lattice. 



an infinitely degenerate ground state at H < H sgut is ac- 
companied by condensation of a macroscopic number of 
zero-energy modes.ollj This should be contrasted with 
the behavior of non- or weakly-frustrated magnets, where 
the antiferromagnetic state below iJ sa t is described by a 
certain wave-vector and, therefore, the phase transition 
at H = H sll t corresponds to condensation of one or a few 
zero-energy modes. In the vicinity of -ff S at the thermo- 
dynamic properties of a strongly frustrated magnet are, 
therefore, controlled by flat branches of excitations with 
vanishing energy at H = H S3 ±- 

The following theoretical analysis is divided into two 
parts. The next section gives a general thermodynamic 
description of magnets with a macroscopic number of 
soft excitations. We also calculate the number of such 
modes for the three frustrated spin models. The sub- 
sequent section describes Monte Carlo results for vari- 
ous thermodynamic quantities of frustrated models and 
also simulations of the adiabatic demagnetization process 
for a few frustrated and non-frustrated magnets. These 
are followed by discussion of suitable magnetic materials, 
which can exhibit the predicted behavior. 



II. THEORY 

Our aim in this section is to understand the effect of 
soft modes on a cooling rate of a strongly frustrated mag- 
net under adiabatic demagnetization. The cooling rate 
is related to isothermal characteristics by 



dT 
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(2) 



where C is the specific heat in a constant magnetic field. 
We shall also use the normalized rate dividing (|J) by 
the cooling rate of an ideal paramagnet (dT / dH) p ^ ra = 
T/H. Values (dT/dH)^ mm > 1 indicate that a magnet 
cools down faster than a paramagnet in a certain range 
of fields and temperatures. The cooling rate is enhanced 
for systems with large negative values of (dS/8H)T- 



A. Magnetothermodynamics 

A classical magnet on an JV-site lattice has 2N degrees 
of freedom in total with two modes per spin. In the sat- 
urated phase above H sat all 2N magnetic excitations are 
gapped. In the harmonic approximation the excitation 
energies of non-frustrated magnets are given by 



w(k) =H- ff sat + e(k) 



(3) 



where e(k) is nonnegative and vanishes at one or a few 
points in the Brillouin zone. Non-Bravais lattices have 
several modes for every wave- vector k and below we im- 
ply summation over all such branches. At low tempera- 
tures anharmonic effects are small, which allows to find 
explicitly the fluctuation contributions to various phys- 
ical quantities. In particular, the specific heat per one 
spin is C — 1 (fee = 1) and 
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Since e(k) ~ k 2 near its zeroes, (dS/dH)T at H = H sat 
is a T-independent constant in 3D and has a weak loga- 
rithmic singularity in 2D. 

The excitation spectra of strongly frustrated magnets 
in the saturated phase are special in that e(k) = in 
the whole Brillouin zone for several flat branches of ex- 
citations. We denote the number of such modes as Na, 
N2 being the number of ordinary dispersive modes (|^) 
and N2 + N4, = 2N. As field approaches H sat , the con- 
tribution from soft dispersionless excitations grows and 
the harmonic approximation fails. To study the effect of 
soft modes one needs to consider an anharmonic energy 
functional of the type: 

E = J2(H - H sat )\y k \ 2 + ]T V£|£ +q y;_ q tf P Wc, (5) 

k k.p.q 



where y^'s denote collective coordinates of soft modes 
and V^ p ~ J are nonnegative interaction constants. Cu- 
bic terms vanish in the expansion by a spin-rotational 
symmetry. Also, an interaction between soft and the har- 
monic modes can be neglected at low temperatures. The 
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above energy functional determines scaling of the corre- 
sponding contribution to the free energy: 



F 4 = -^T\nT-Tf 

4 J \VT 



(6) 



The function f(u) has the following asymptotes: f(u) as 
/(0) + uf'(0) with /(0) > and /'(0) < as u -> 0, and 
/(u) ~ —(N4/2) In it as it — ► 00. 

It follows from Eq. (^) that the specific heat and the 
magnetic susceptibility depend on magnetic field—and 
temperature only via combination (H — if sa t)/vT'0 In 
particular, the specific heat per one spin reaches a univer- 
sal fractional value at H = H sat : C = (^N 2 + jN 4 )/N = 
l — N4,/(AN). For the magnetocaloric effect, Eq. (||) gives 
at H = H sat : 
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dH 



N 4 



(7) 



which significantly exceeds the contribution of dispersive 
modes Eq. Thus, condensation of a macroscopic 

number of zero-energy modes at the saturation field 
greatly enhances the magnetocaloric effect in strongly 
frustrated magnets at low temperatures. An important 
characteristics of a magnet is, then, the number of soft 
modes N 4 . 



B. Soft modes 

1. Kagome antiferromagnet 

In order to find the number of soft modes we perform 
classical spin-wave calculations in the saturated collinear 
state. A primitive unit cell of a kagome lattice contains 
one triangular plaquette. The corresponding Bravais lat- 
tice has N/3 states. Expanding spin components along 
the field in small transverse deviations we obtain a clas- 
sical spin-wave Hamiltonian to a desired order. In the 
harmonic approximation two transverse polarizations for 
the spin excitations are independent. As a result, the 
spectrum consists of three double-degenerate modes: 



wi(k) = H-6J , 
W2 )3 (k) = H-3J± J y /3{l + 2 lk ) 



(8) 



where 7k = ^(cosk x + 2 cos \k x cos ^k y ). The total 
number of soft quartic modes at H saii — 6J is, there- 
fore, N 4 as = 2N/3. This is exactly twice the number 
of hexagon voids on a kagome lattice. The correspon- 
dence between such voids and soft modes has firsUbeen 
found for a kagome antiferromagnet in zero_fielalj and 
was, later, extended to finite magnetic fieldsO Local soft 
modes in the collinear state at H — H sat consist in al- 
ternate tilting of spins around elementary hexagons in 
two different directions perpendicular to the field. The 



geometric picture of local soft modes can be straightfor- 
wardly generalized for a classical spin model on a checker- 
board lattice, which is the other planat-frustrated lat- 
tice with voids around empty squares. cB The checker- 
board antiferromagnet has N local zero-energy modes at 
Hsat = 8 J. 



2. Pyrochlore antiferromagnet 

The spin-wave calculations for a pyrochlore antiferro- 
magnet yield four double-degenerate modes defined in 
the Brillouin zone of a face-centered cubic (fee) lattice, 
which has N/4 states: 



wi f2 (k) = H-8J , 

w 3 , 4 (k) = H - 4J ± 2 jyi + 3^ 



(9) 



i (cos k x + cos k y + 4 cos k z 



COS 2 
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The total number of zero modes in a pyrochlore antiferro- 
magnet at the saturation field is N^ yi = N, i.e. the same 




(a) 




(b) 



FIG. 2: Two projections of a pyrochlore lattice on a (001) 
plane. Small, light and dashed, squares represent one type of 
tetrahedra; big squares with crosses correspond to the other 
type of tetrahedra. Small light tetrahedra belong to the same 
(001) plane of an fee lattice. Small dashed tetrahedra lie in 
an adjacent plane: (a) below and (b) above. Soft spin modes 
on hexagon loops Z123 and £134 are shown with +'s and — 's on 
the top panel. 
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as for the checker-board model. The geometric interpre- 
tation of these modes becomes, however, more compli- 
cated. A dispersionless soft mode can still be attributed 
to a local perturbation in the collinear structure: alter- 
nate tilting of spins around a closed line, which shares 
two vertices with each crossed tetrahedron. The shortest 
such loop on a pyrochlore lattice is a hexagon. Simple 
counting yields N hexagons for an iV-site pyrochlore lat- 
tice. If a one-to-one correspondence between the soft 
modes and the hexagons holds, then the above number 
would imply 2N soft modes, which is twice more than 
the spin-wave prediction. The contradiction is resolved 
in favor of the spin- wave result by establishing that some 
of the hexagon loops on a pyrochlore lattice are linearly 
dependent with total N/2 linear relations. This brings 
the number of independent hexagon spin-loops (modes) 
to N/2 in agreement with the spin-wave expansion. Be- 
low, we briefly sketch a proof of the above statement. 

A pyrochlore lattice, Fig. [l]c, is formed with two types 
of tetrahedra. To simplify geometric representation we 
draw a pyrochlore structure as an fee lattice constructed 
from one type of tetrahedra. A cut through such a lattice 
along (001) plane is shown in Fig. ^. Hexagons in Fig. || 
correspond to closed loops which include two tetrahedra 
in one layer and one tetrahedron in an adjacent layer. 
Soft spin modes on two such loops are shown in Fig. ^ 
with pluses and minuses, which indicate alternate tilting 
of spins from the field direction. The first type of a linear 
relation between spin-loops appears for a block formed by 
four tetrahedra numbered 1-4 in Fig. |^a. A superposition 
of the local soft modes corresponding to hexagon loops 
Zi23 and /134 leaves two inner spins unmoved and tilts 
only outer spins around the block perimeter. The same 
eight-spin mode is obtained by a superposition of spin- 



loops Z124 and ^234- Thus, every such block of four tetra- 
hedra provides one linear relation between four loops. 
This reduces the number of independent spin-loops by 
the number of such blocks iV/4. The second type of a 
linear relation corresponds to a block of six tetrahedra: 
four in one layer, 1, 3, 5, and 6, and two in adjacent lay- 
ers, 2 and 7, Figs. |]a,b. A superposition of spin-loops 
Z123 and ^256 corresponds to the same mode as a super- 
position of Z157 and l^j. Counting all such blocks brings 
additional A/4 linear relations. This makes in total N/2 
linear relations between hexagonal spin-loops and, hence, 
A| yr = 2{N - N/2) = N soft modes for a pyrochlore 
antifcrromagnet at the saturation field. Composite low- 
energy excitations on six-spin loops, which are similar 
to the above soft modes, have been recently observed in 
neutron scattering studies of pyrochlore antiferromagnet 
ZnCr 2 4 in zero fieldO 



3. Garnet antiferromagnet 

The structure of a garnet lattice is by far the most com- 
plicated among frustrated lattices. In GGG the magnetic 
Gd 3+ ions form two interpenetrating triangular sublat- 
tices with 24 atoms per a cubic unit cell. Exchange inter- 
action couples each spin with its nearest-neighbors from 
the same sublattice. Spins on two different sublattices 
interact only via weaker dipolar forces. In the present 
study we focus on the behavior of an exchange model and, 
therefore, consider only one sublattice of corner-sharing 
triangles with 12 spins per a unit cell represented, e.g., by 
dark circles in Fig. |l|b. The nearest-neighbor exchange 
Hamiltonian on a garnet lattice can be written as 



W — J [Sij • (S 2 i + S 3i + S ii+X + S 5i+X ) + S 2 i ■ (S 3i + S 6i+y + S 7i+y ) + S 3i • (S 8i+2 + Sg i+2 ) 

i 

+ &4i ■ (Sgj + S$i+z + Sio;+ 2 ) + Ssi ■ (Sgi-a: + Sih) + S(a ■ [Sn + Snj+a) + S 7 i ■ (Sgi-y + Sl2j) 



Ssi ■ (Sgi + Sioi) + Sg; • Si2i+J/ + Sioi • (Slij + Si2j) + Sn, • Si2i 



Hz 



ccman ; 



(10) 



where the first spin index j = 1-12 is the site number in 
a unit cell and the second spin index i is the cell number 
on a cubic lattice; i±x,... being adjacent cells along three 
orthogonal directions. 

After Fourier transformation the spin-wave energies 
are given by eigenvalues of a 12 x 12 matrix. The cor- 
responding equation factorizes into a product oiJ&j — 
H + 6J) 4 and an eighth-order polynomial in w£3 At 
H = 6J there are exactly four zero roots for an arbi- 
trary value of the wave-vector. Zero eigenvalues deter- 
mine 7Vf ar = 2N/3 soft modes at the saturation field, 
which is the same number as for a kagome antiferromag- 



net. Thus, the term Ziyper-kagome used sometimes to 
refer to a garnet structure receives an extra justification. 
Geometric picture of soft modes in a hyper-kagome an- 
tiferromagnet is, however, much more complicated than 
in its planar counterpart. The shortest closed loop on a 
garnet lattice, which corresponds to an elementary soft 
mode, crosses ten triangles. At the same time each bond 
between a pair pf nearest-neighbor spins participates in 
five such loops.tj This should be contrasted with sim- 
ple hexagon loops on a kagome lattice, where every bond 
is included in one loop only. The total number of ten- 
triangle loops for a garnet lattice is calculated to be N. 
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FIG. 3: Scaling behavior of the specific heat in the satu- 
rated phase for a kagome (top panel) and a pyrochlore (bot- 
tom panel) antiferromagnet. Different symbols correspond to 
field-scans performed at several fixed temperatures. Arrows 
indicate the limiting values of the specific heat per one site 
at the saturation field. Fields and temperatures are given in 
units of J. 



A comparison with iV/3 soft modes for one transverse po- 
larization found in the spin-wave analysis suggests that 
there are 2N/3 linear relations between closed loops for 
a garnet structure. An underlying geometric picture be- 
hind these relations remains to be clarified. 

To conclude this section we note that a pyrochlore 
lattice is the most frustrated one among the above 
three structures, because it has the largest number 
of soft modes. The same is also true in zero mag- 
netic field, where a classical Heisenberg pyrochlore lat- 
tice antiferromagriet remains disordered down to low- 
est temperatures,E2l whereas a classical kagome an- 
tiferromagnet— .develops a unique planar triatic order 
parameter Mi 1 ^ As a result, a pyrochlore antiferromag- 
net should exhibit the fastest cooling rate under adiabatic 
demagnetization among frustrated magnets. 



III. MONTE CARLO SIMULATIONS 

Monte Carlo (MC) simulations of the three spin mod- 
els have been performed with the standard Metropolis 
algorithm. To improve statistics at low temperatures we 
introduce a maximum change AS Z = 3T for a choice 
of a new spin direction in the local coordinate frame. 
In this way an acceptance rate for spin moves is kept 
at the level of 30%. (Other values AS Z /T ~ 1-5 will 
work as well.) MC simulations are done on finite clus- 
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FIG. 4: Scaling behavior of (dS / 8H)t at H = H aa ,t in kagome 
and pyrochlore antiferromagnets. Inset shows the field depen- 
dence of (dS I 8H)t at T = 0.01J for a pyrochlore antiferro- 
magnet. Fields and temperatures are given in units of J. 



ters with periodic boundary conditions. Typical clus- 
ter sizes are N = 3 x 18 2 = 972 for a kagome lat- 
tice, N = 4 x 6 3 = 864 for a pyrochlore lattice, and 
N = 12 x 5 3 = 1500 for a garnet lattice. Compari- 
son with smaller and larger clusters have been made to 
insure that finite-size corrections are small. Typical MC 
runs are carried with 10 5 MC steps per spin for equilibra- 
tion followed by averaging over 10 5 measurements made 
in intervals of 5-10 MC steps. Statistical errors are es- 
timated by performing several (~ 5) runs from different 
initial configurations. The errors for all presented data 
are small and do not exceed symbol sizes. The specific 
heat is calculated from the variance of the total energy: 
C = ((E 2 ) - (E) 2 )/T 2 , while the magnetocaloric effect 
is expressed via a cumulant of the energy and the mag- 
netization: (dS/dH) T = ((EM) - (E)(M))/T. 

We present in Fig. || the specific heat per one site of 
kagome and pyrochlore antiferromagnets plotted with re- 
spect to the scaling parameter u = (H — H sa t)/VT > 0. 
The data have been collected during field-scans at fixed 
temperatures. Points from different scans, except of the 
highest temperature of T = 0.1J, fall on the same curve 
confirming the predicted scaling form. The specific heat 
at H = iJ Ba t (u = 0) reaches values which are close to 
the asymptotic results derived from the number of soft 
quartic modes: C = 5/6 and C = 3/4 for a kagome and 
a pyrochlore antiferromagnet, respectively. The magne- 
tocaloric effect near saturation is presented in Fig. 0. The 
main plot shows temperature dependence of (dS/dH)T 
exactly at the saturation field for kagome, pyrochlore and 
square lattice (576-site cluster) antiferromagnets. Nu- 
merical results for a garnet antiferromagnet are not in- 
cluded because at H > H sgut they virtually coincide with 
the data for a kagome antiferromagnet. Significant dif- 
ference (Fig. ||) between the frustrated magnets on one 
side and a non-frustrated square lattice antiferromagnet 
on the other side confirms an enhancement of the magne- 
tocaloric effect related to magnetic frustration. A range 
of magnetic fields in the vicinity of H sat , where such an 
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FIG. 5: Cooling rates of pyrochlore and square lattice anti- 
ferromagnets normalized to a cooling rate of an ideal param- 
agnet. Fields and temperatures are given in units of J. 



enhancement occurs, is illustrated on the inset. 

Direct comparison of the cooling rates of magnets with 
and without frustration is made in Fig. [^, where results 
for a pyrochlore and a square lattice antiferromagnet are 
plotted for several field scans. The two magnets have the 
same saturation field i/ sa t = 8J. For a pyrochlore an- 
tiferromagnet this field corresponds to a crossover from 
the nondegenerate saturated phase into a highly degen- 
erate state below i? sa t . The crossover is accompanied by 
a large increase in the adiabatic temperature variation, 
which exceeds the cooling rate of an ideal paramagnet 
by more than an order of magnitude. In contrast the 
spin model on a square lattice has a real transition near 
-ffsat into a nondegenerate antiferromagnetic state with 
divergent correlation length. The cooling rate becomes 
negative below the transition, which means that a mag- 
net heats up during demagnetization rather than cooling 
down. There is also no significant enhancement of the 
cooling rate above -ff S at apart from a small narrow peak 
due to an anomaly in (dS '/ 8H)t at the transition. 

Finally, a full adiabatic demagnetization process can 
be simulated with MC by measuring at fixed T and H 
isothermal characteristics C and (8S/dH)T to find the 
cooling rate (^|) and, then, choosing a new value of tem- 
perature according to AT a d/AH = (dT/dH)s for suf- 
ficiently small AH. Adiabatic temperature changes for 
demagnetization from two starting temperatures T = J 
and T = 0.1J and the same initial field H — 9J are shown 
for garnet, pyrochlore, and cubic (N = 512) antiferro- 
magnets in Fig. [|. Since saturation fields for the three 
magnets are different, we use rescaled exchange constants 
J* = 0.75 J for a pyrochlore and J* = 0.5 J for a cubic 
antiferromagnet to set the common saturation field of 
-ffsat = 6J in all three cases. During demagnetization 
the non-frustrated cubic antiferromagnet has the small- 
est relative temperature variation and starts to heat up 
below the transition into antiferromagnetically ordered 
state. The adiabatic temperature decrease for frustrated 
magnets is much larger with the largest effect exhib- 
ited by a pyrochlore antiferromagnet. Note, that tem- 



FIG. 6: Temperature variation during adiabatic demagneti- 
zation of a pyrochlore (full line) , a garnet (dashed line) , and a 
cubic (dot-dashed line) antiferromagnet. Dotted lines repre- 
sent demagnetization of an ideal paramagnet (T/H — const). 
Fields and temperatures are given in units of J. 

perature of a pyrochlore antiferromagnet drops by more 
than ten times as field decreases from H[ n [ = 1.57J sa t to 
Hfi n — 0.8-f/sat- To achieve a similar cooling with a para- 
magnetic salt magnetic field has to be decreased at least 
by a factor of ten. The two distinctive features of the de- 
magnetization of a pyrochlore antiferromagnet are (i) a 
persistent temperature decrease up to a zero field, though 
with a reduced cooling rate below -ff sa t) and (ii) an in- 
crease in the relative temperature variation as a starting 
temperature goes down, which is a consequence of the 
low-temperature singularity (0). 



IV. DISCUSSION 

The practical limitation on experimental observation of 
an enhanced magnetocaloric effect comes from a strength 
of the saturation field. At present, only frustrated mag- 
nets with exchange constants of the order of few Kelvin 
or less can be saturated with dc magnetic fields. The 
primejacample is gadolinium gallium garnet with H sat ~ 
1.7 T,i3 which is a popular refrigerant material in an jap* 
plied experimental research on the magnetic cooling. clE-II 
Quantitative comparison of the cooling rates of an ideal 
paramagnet and a Heisenberg garnet antiferromagnet is 
made in Fig. ^. Starting with the same temperature a 
garnet antiferromagnet cools down to much lower tem- 
peratures than a paramagnet in a wide range of mag- 
netic fields. There are also two py-rqch lore antiferromag- 
nets Gd 2 Ti 2 7 and Gd 2 Sn 2 7 |aH3'El which have satu- 
ration fields about 7 T. It is an interesting experimental 
challenge to observe the fast temperature decrease upon 
adiabatic demagnetization in these materials. Operation 
with relatively high magnetic fields makes rather ques- 
tionable a technological application of the two gadolin- 
ium pyrochlores. The theoretical results suggest, though, 
a new direction for a search of suitable materials: a py- 
rochlore antiferromagnet with a factor of four smaller 
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exchange constant than in Ga^T^Oy would be a more 
efficient refrigerant than Gd3Ga 5 0i2, while operating in 
a similar range of fields and temperatures. 

Real materials have, of course, more complicated mag- 
netic interactions than just a nearest-neighbor exchange 
between classical spins. For example, in GGG a dipole- 
dipole interaction between neighboring spins anmunts to 
20-30% of the corresponding exchange energy£3 Long- 
range dipolar interactions play «1sp a significant role in 
the two pyrochlore compounds .Ej€3 Influence of these 
and other extra perturbations on the magnetocaloric ef- 
fect in frustrated magnets requires further theoretical 
analysis. Meantime, one can make several general re- 
marks on the relative importance of various perturba- 
tions. The quantum effects omitted in the present study 
do not destroy the anomalous magnetothermodynamics 
near H sat - Similar to their classical counterparts quan- 
tum antiferromagnets on the three-frustrated lattices also 
have flat branches of excitations.Ej The anomaly (0) is, 



nevertheless, affected by a weak dispersion (~ Jdisp) of 
soft excitations induced, e.g., by a second-neighbor ex- 
change or dipole-dipole interactions. The magnetocaloric 
effect will saturate below T* ~ Jdisp at (dS/dH)r ~ 
1/vT*) which could be relatively large. Such a frus- 
trated magnet can still be used, though less efficiently, 
as a refrigerant to reach temperatures below T*. A fur- 
ther analysis of a destructive role of various perturbations 
should focus, therefore, on an induced dispersion in flat 
magnon branches. 
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